function coor = findbrightest(I)
    [sx,sy,~] = size(I);
    h1 = fspecial('gaussian',[sx sy],(sx+sy)/10);
    hsv = rgb2hsv(I);
    v = hsv(:,:,3);
    h = fspecial('disk',4);
    v = imfilter(v,h);
    g = gradient(double(v));
    g = g .* h1;
    max_g = max(g(:));
    [r c] = find(g == max_g);
    coor = [c r];
end